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I. INTRODUCTION 

Let us consider the Laplace operator on functions over some domain M C IR" with 
Cartesian coordinates x°' — , . . . , x"^: 

A = g'^'dad,, 5a = (1.1) 

with a constant positive definite metric (gab) = ig^^)~^, and the matrix Schrodinger oper- 
ator 

L = IA + U, (L2) 

where / is the N x N unit matrix and U = U{x) some N x N matrix- valued potential. We 

assume the potential U to be Hermitian, = U, and smooth in the domain M. If M is 
noncompact, it is assumed, in addition, that U{x) decreases sufficiently rapidly at infinity. 

It is well known that L is an elliptic and self-adjoint operator with a real spectrum 
bounded from above on a Hilbert space of functions over M with some suitable boundary 
conditions. In the case of a periodic potential the boundary conditions can be replaced 
by periodicity conditions, and in the case of a noncompact domain they take the form of 
some asymptotic conditions at infinity, requiring sufficiently fast decrease. 

This implies that for any t > the operator K(t) — exp(tL) is well defined. Its kernel 
K{x,x';t), often referred to as the heat kernel, is the fundamental solution of the heat 
equation 

{Idt - L)K = 0; (1.3) 
it satisfies a singular initial condition 

K{x, x'; +0) = {detgab)-^^^S{x - x')I, (1.4) 

where 5{x) denotes the Dirac delta-distribution. It is well known that when the points x 
and x' do not lie on the boundary, there holds an asymptotic expansion as t — > -|-0 near 
the diagonal x x' oi the form 

K{x,x';t) ~ {47rt)-''/^expi-^\x-x'\AY2t''Hkix,x'), (1.5) 

^ fe=0 

where |a;| = gabX°'x^- Prom (1.3) and (1.5) it follows that the sequence of the coefficients 
Hi~ is characterized by the differential-recursion relations 

{D + k) Hk = LHk-i, ioTk>l, (1.6) 

where we abbreviate 

D={x-xTda, (1.7) 

together with the initial condition 

Ho = I. (1.8) 
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The analysis of the heat kernel, or Hadamard, coefficients = Hk{x, x'), (k = 
0, 1, 2 . . .), introduced by (1.6)-(1.8) is the subject of the present paper. Their role goes far 
beyond the heat kernel. They can be introduced for any linear second order differential 
operator with a Laplace-like principal part (cf. e.g.^~^). The sequence {Hk) turns out to 
have diverse mathematical and physical applications: spectral geometry, Hadamard's ele- 
mentary solution, Huygens' principle, statistical physics, determinants of elliptic operators, 
anomalies, Korteweg-de Vries hierarchy of soliton equations etc. 

In the present paper we assume the metric (gab) to be Euclidean, i.e. positive defi- 
nite, and the operator L to be elliptic for the sake of simplicity. The formal results hold 
for a pseudo-Euclidean metric of any signature and hyperbolic (normally hyperbolic or 
ultrahyperbolic) operators too. 

Our paper is organized as follows. In Sects. II and III we obtain explicit formulas for 
the Taylor coefficients [d^Hklip = 0,1,2,...) of Hk{x,x'), where x is the variable point 
and x' is the origin. We analyse these quantities and determine some special terms and a 
'common denominator' (see below). In Sect. IV we consider the one-dimensional case and 
derive a simplified algorithm for the diagonal values [Hk] of the via a matrix Riccati 
equation. Then we generalize the so-called 'Lenard recursion' from the scalar case to the 
matrix one. The new recursion is formally non-local, but can be effectively done in the 
space of differential polynomials in U. It leads to a proposal of some matrix Korteweg-de 
Vries hierarchy for U = U{x,t). 

It is not possible to quote all the literature on the heat kernel approach. Let us 
mention papers which deal with the matrix case^~^^. The methods applied there differ 
from ours. 

II. THE HEAT KERNEL COEFFICIENTS 

For a smooth potential U = U{x) the coefficients Hk{x,x') are known to be smooth 
near the diagonal x = x'. Thus, to study the sequence (-ff/c) we use the Taylor expansion 
with respect to the variable x at the origin x' taken as an asymptotic series as a; — > x': 

oo ^ 

Hk{x,x') ~ V -(x - x'r^ ■■■{x- x'rnda, ■ ■■da.Hk]. (2.1) 

Here and in the following the square brackets [ ] denote the diagonal value of a two-point 
quantity H = H{x, x'): 

[H] = [H]{x) := H{x,x). (2.2) 

The coefficients of the Taylor expansion can be arranged to matrix-valued symmetric 
differential forms (or covariant tensors), shortly p-forms,. The symmetric product V of two 
such forms A^p^ and B^^^ is defined as 

V = A^ai-apBb^-\)dx°'^ (8) • • • (8) dx"-^ (g) dx^"^ (g) • • • (g) dx^", (2.3) 

where the parentheses mean symmetrization over all indices included. The symmetric 
differential 

d : ^(p) ^ ^(p+i) (2.4) 
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on these forms is defined by 

= d^a^Aa^.^a^^^^dx"-^ dx'^^+K (2.5) 

Using these notations, we have 

[d^Hk] = [da, ■ ■ ■ da^Hk]dx''' • • • (8) dx''-. (2.6) 
Further, we define some trace operation 

g-^ : ^(p) ^ A(p_2) (2.7) 

by 

g-^A^o) = 0, ^"'^(1) = 0, 

g-^A^p^ = g'''^'A^a,a,as...a,)dx''' ® • • • ® dx''- , for p > 2. (2.8) 

The Laplacian on symmetric forms 

A := g-^d^ : A^^^ ^ A^^^ (2.9) 

generalizes the usual Laplacian (1.1) on functions, i.e. 0— forms. 

Let us now apply d^ to (1.6) and then restrict to the diagonal x = x'. Taking into 
account the commutation relation 

[dP,D]^pdP, (2.10) 
we obtain the purely algebraic recursion system 

(p + k)[dPHk] = g-^[dP+^Hk-i] + (^) Uip-q) V [d'^Hk-i], for /c > 1, p>0 

(2.11) 

where f/(p) — d^U, with the initial conditions 

[Ho] = /, [dPHo] =0 for p > 1. (2.12) 

An explicit non-recursive solution of (2.11), (2.12) can be presented by means of the 
following symbolism. Let us define linear operators Mq, [q = —2, —1, 0, 1, . . .) on matrix- 
valued symmetric forms according to 

M_2^(p) := fir-'^(p), M_i^(p) := 0, (2.13) 
M,yl(p) + V A(p), for q > 0. (2.14) 
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Then the system (2.11) can be rewritten in a short-hand notation: 

[dPHk] = -—TMp-,[d'^Hk-i]. (2.15) 

By applying mathematical induction to this recursion system we obtain the following. 
Theorem 2.1. There holds for k > l,p > 

92, \r=l ' / 

where the sum runs over the integers q2,Q3, ■ ■ • ,Qk so that 

< Qr+i < ?r + 2, for r = 1, 2, . . . , /c — 1, 

and 

This result can also be obtained, as a special case, from the general explicit formula for 
the heat kernel coefficients for any second-order elliptic operator of Laplace type obtained 
by one of the authors (I.G.A.) in^. 

The formula (2.16) is not well suited to recognize special polynomial constituents of 
[d^iiffc]. To do this, we decompose [d^H]^] into its homogeneous parts of degree d with 
respect to U : 

[d^iy,] = X:<]. (2.17) 

The system (2.11) decomposes in an easy way into homogeneous parts too. Thereby 
the linear and the highest degree parts decouple from the whole system: 

iP + k)H^^l-9-'H^X% (2-18) 
iP + k)Hii - E h U^,-,) V Hi\,_,. (2.19) 

Q=0 



The initial values are given by 

, -,^rr(p) r, 

'{P) 

Mathematical induction gives the following. 
Theorem 2.2. The hnear part of [d^H]^] equals 



{p+l)HS^U^py (2.20) 



^i^''- {p+2k-i)r ^ ^' ^^-^^^ 
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while the highest degree part of ujfl of [d^Hk] is 

where the sum runs over the integers q2,Q3, ■ ■ - Qk, so that 

< Qr+i < Qr for r = 1,2, . . . , k — 1 

Qi = P, Qk+i = 0. 

Let us call a number-theoretical function d{p, k) a comm.on denominator to the system of 
differential polynomials [d^H]^ if all the modified polynomials in U, U(^i^, U{2)i • • ■ 

H^^^ :=d{p,k)[dPHk] (2.23) 

have integer coefficients. 

Proposition 2.3. A common denominator in the aforesaid sense is given by 



d{p, k) = -M for even k = 21, (2.24) 



and 



d{p, k) = ^U + for odd k^2l-l, (2.25) 



j=0 

Proof. Substituting 



into (2.11) we get the recursion system for the Hjf'' 

- (p) _ d{p,k) ^-1 fr(P+2) ^ V" M d{p,k) -(g) 



{p + k)d{p + 2,k-iy ^W(p + '^K?>^-i)^''"'^ 



fc-i' 



(2.27) 

for A; > l,p > 0. One can show that the coefficients of this system are integer, which 
proves the theorem. 

To find the least common denominator remains an open problem. Fulling^^ found the 
expressions (2.24), (2.25) except for the factor l/p\. It is obvious, that this factor does not 
matter for p = 0. 
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III. THE DIFFERENTIAL POLYNOMIALS [Hk] 

Let us analyse more closely the structure of the differential polynomials [Hj^]. The 
common denominator for the system [H^], {k — 1,2, ■ ■ •) reads 

n '^^li^- (3-1) 

The extremal terms are given by 

where the points indicate terms of a degree higher than 1 and less than k. The sequence 
begins with 

[^^i] = u, 

[H2] = l{AU + 3U^), (3.3) 
[Hs] = {A^U + 5 [U{AU) + {AU)U + g-\dU V dU)] + 10C/=^} . 

Using the identity 

2g-^{dWi V dW2) = A{WiW2) - {AWi)W2 - Wi{AW2), (3.4) 
whith Wi, W2 being 0-forms , one can eliminate the operator g~^ here: 

[H^] = ^ [2A^U + W{AU) + h{AU)U + 5A(t/2) + 20t/^} . (3.5) 

J. 

Let us call a monomial or differential polynomial simple if it is solely composed from 
U and A (but not from d and g~^), and the total number of symbols U and A in a simple 
monomial the weight w. It is obvious that a simple monomial is always a 0-form and that 
the number of symbols U in it equals its degree. 

Theorem 3.1. There are exactly 

simple monomials of degree d and weight w. 
Proof. A simple monomial W reads W = or 

W = U'^°{AWi)W2, (3.7) 
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where ?7^° and W2 may be empty, while Wi is not. Since the decomposition is unique, the 
number C{w,d) of monomials W follows from the numbers C{wi^ di), C{w2i c/2) of W^i, W2 
respectively. Thus, 

C{w,d)^Y^ C{wi,di)C{w2,d2) (3.8) 

do=0 di,d2;wi,W2 

for d < w, where the inner sum runs over integers wi,W2, di, d2 so that 

wi > di > 1, W2 > d2 > 0, wi + W2 = w — do — 1, di + d2 = d — do. 
The recursion formula (3.8) can be translated into a quadratic equation 

xyF^ + {x + y-l)F+l = (3.9) 

for the generating function 

00 w 

F = Fix, y) = Y.Yl '^(«'' '^)^'"'^™~^ (3-10) 

w=l d=l 

with the initial condition 

F(0,0) = 1. (3.11) 

Further, using 

(i+^r-iE^C:)^™ (3.12) 

w=0 ^ ' 

we expand the solution of the equation (3.9) 

F{x, y) = l(\-x-y^ ^{X-x-yf-^x^ (3.13) 

in the power series (3.10) and prove that (3.6) is correct. 

B. Fiedler^'^ derived special cases of theorem 3.1, namely for w = 2 and = 3, by 
combinatorial arguments. 

Now, let us try to eliminate from the differential polynomials by means of the 
identities like (3.4). Making use of 

2g-^{dWx V V dW^) = A{WiW2W3) - {A{WiW2))W3 - WiA{W2W3) + Wi{AW2)W3. 

(3.14) 

where Wi are 0-forms, we easily get a simple form for [114] : 

[H4] = I ^A^U + ^ [A'^U^ + {A^U)U + UA^U + {AUf + A{{AU)U) + A{U{AU))\ 
+ A(C/3) + (A(C/2))C/ + C/A(C/2) + (AC/)C/2 + U'^AU + U{AU)U + 5U^\. (3.15) 
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However, it is impossible to do this for higher polynomials [Hk],k > 5, because of terms 
like 

One of the authors (R.Sch.) has calculated [Hi], . . . , [H^] in simple form for the matrix 
Schrddinger operator in the early paper^^. Later T. Osborn et al.^ obtained the result in 
a "non-simple" form. The other author (I.G.A.) calculated the coefficients [Hi], . . . , [H4] 
in the general case, that means for every Laplace-type operator, in^^'-*^^'^. 



lY. THE ONE-DIMENSIONAL CASE 

Let us discuss more detailed the one-dimensional Schrodinger operator 

L = cf + U, (4.1) 
where now d = d/dx and n=l. In this case the recursion system (2.11), (2.12) reduces to 



{p + k)[dPHk] = [d^+^Hk.i] + J2 (^) Up-S'Hk- 



for/c>l,j9>0 (4.2) 



g=0 

[Ho]=I, [d^Ho]=Q, forp>l. 

The explicit solution (2.16 ) of this system can be rewritten now in a better form. 
Theorem 4.1. There holds for k >2,p>0 

qi,---,qk r=i ^ 

where 

t/_2 := /, := 0, Up := Up = dW for p > 0, (4.4) 

c{p,q):^{^^+5l+\ (4.5) 

and the sum runs over the integers q2,Q3, • • • ,Qk so that 

< Qr+i < Q'r + 2 for r = 1, 2, • • • , /c — 1, 

Qi = P, Qk+i = 0. 

A proof has been given by R.Sch. in^^ for the scalar case (A^ = 1). Inspection shows 
that the arguments of^^ do not depend on the number A^ of components and are valid for 
the matrix case too. FuUing^^ found that formula (4.3) is well suited for calculating the 
quantities [d^Hk] by means of computer algebra. This theorem also follows, as a rather 
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special case, from the results obtained by I.G.A. in^, where a general explicit formula for 
the heat kernel coefficients to any second-order elliptic operator of Laplace type is obtained. 

From Theorem 4.1. it follows that the quantity [d^Hi^] for k > l,p > can be brought 
into the form 

k 

[dPHk]=J2 Yl a{ri,r2,...,rd)Ur,Ur,---Ur„ (4.6) 

d=l ri,...,rd>0 

where a(ri, r2, . . . , r^) are some positive coeflBcients and the inner sum runs over the inte- 
gers ri,r2 ■ ■ ■ ,rd so that 

ri+r2-\ h rd = p + 2{k - d). 

There is a one-to-one correspondence between the ordered partitions (ri + 1, r2 + 
1, ■ ■ ■ , + 1) of s — X^^^i(r'i + 1) — p + 2k — d into d numbers and the monomials 
Ur^Ur2 ■ ■ ■ Ur^- Siucc there are exactly (^l]^) such partitions, the total number of monomi- 
als, i.e. the total number of the terms in the sum (4.6) equals 



k 



d=l 



d-l 



The formula (4.3) is explicit, i.e. non-recursive, but somewhat complicated because 
each coefficient c(gr-) 9r-i-i) comprises the two cases qr+i < and qr+i = Qr + 2. As an 
alternative, we offer a recursion for the sequence [Hk], {k = 0, 1, 2, . . .) which is by far more 
effective than the double recursion (4.2) for the double sequence [d'^Hk], {p, /c = 0, 1, 2, . . .). 

Theorem 4.2. Let the differential polynomials Z^, {k = 1,2, . . .) in U = U{x) be recur- 
sively defined by 

k-l 

Zk+l = dZk + ^ ZmZk-m (4.8) 
m=l 



Then tlie quantities 



2k\ 

are determined by tlie recursion 



Zi = U. (4.9) 



Gk := ^[Hk] (4.10) 



Gfe = 2 ^ Z2m-lGk-m: for k > 1, (4.11) 
m=l 

together with the initial condition 

Go = \l. (4.12) 
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Proof. Let us start with the heat kernel K{x, x'; t) to L. Its Laplace transform 

oo 

G{x,x';X)= J dtexp{-tX)K{x,x';t). (4.13) 



is the resolvent kernel to L, i.e. the fundamental solution of L — XI, which means 

(d^ + Uix) - XI)Gix, x'; X) = -Six - x')I. (4.14) 

The resolvent kernel is continuous but not differentiable bX x — x' . Near the diagonal, 
X x' , it has the form 

G{x,x';X) = [G{X)]{x') - l\x-x'\ + (x - x')F{x'; X) + 0{{x - x'f), (4.15) 

where [G{X)]{x) := G{x,x;X), so that 

dG{x,x';X) = ~siga{x-x') + F{x';X) + 0{x-x'), (4.16) 

where is the absolute value of x, sign(x) denotes the sign of x and F(x';A) is some 
matrix that does not depend on x. For sufficiently large Re A and x x' the resolvent 
kernel is known to be invertible, so that one can present the derivative of the resolvent in 
the form 

dG{x, x'; A) := Z{x, x'; X)G{x, x'; A). (4.17) 

Then from (4.16), (4.17) we conclude that Z = Z{x,x';X) is discontinuous with a 
jump at a; = a;' that determines the resolvent kernel on the diagonal: 

{Z{x' + 0, x'; X) - Z{x' - 0, x'; A)}[G'(A)] {x') = -I. (4.18) 

Moreover, the equation (4.14) for x ^ x' becomes equivalent to the Riccati equation 
for Z 

dZ{x, x'; A) + Z^{x, x'; A) + U{x) - XI = 0. (4.19) 

Let us look for a solution of this equation in the form of an asymptotic expansion as 
Re A — > oo, |argA| < tt. Using the asymptotic expansion of the heat kernel (1.5) and (4.13) 
we obtain 

Z{x, x'; A) ~ X{x; A) + sign(a; - x')Y{x; A), (4.20) 
where the formal power series 

oo 

X(a;;A) = J](4A)-^Z2fc(x), (4.21) 

k=0 

oo 

Y{x;X) = J](4A)-'=+i/2Z2fc_i(x), (4.22) 



k=0 
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are composed of the quantities Zk{x) defined above. Namely, substituting the expansion 
into the Riccati equation (4.19), we find that Z_i = I, Zq = and Zk{x), (/c = 1,2, . . .), 
are just the diflFerential polynomials in U{x) determined by the recursion system (4.8), 
(4.9). 

Therefore, the asymptotic expansion of Z{x,x'; A) is purely local but discontinuous at 
X = x'. From (4.20) it is clear that the jump of Z{x,x'; A) at a; = x' equals just Y{x; A), 
so that (4.18) takes the form 

2Y{x;X)[G{X)]{x) = I (4.23) 
Using (4.13) and the asymptotic expansion of the heat kernel on the diagonal as t ^ +0 

oo 

[K{t)] ^ (47r)-^/2^t'=-i/2[if,], (4.24) 

fc=0 

we easily get the corresponding asymptotic expansion of the resolvent on the diagonal as 
A — > cxd: 

oo 

[G{X)]^2j2{4X)-'^-'/'Gk, (4.25) 

fc=0 

where 

G.:=^M. (4.26) 

Therefrom, using (4.22) and (4.23), we find that the coefficients G^ are determined by the 
recursion system (4.11), which proves the theorem. 

The result of the Theorem 4.2. can be also obtained from the results of the paper of 
Fulling^, who uses a different method and notation. 

One can explicitly find the linear and quadratic parts of the differential polynomials 
Zk and Gk- 

Theorem 4.3. There holds fork>3 

^fe = f/fe-i + 5^ |(^ ~ j) - l| UpUk-3-p + ■■■, (4.27) 
p=o ^\P'^^ J 

where the points indicate terms of degrees higher than two. 

Proof. The linear i and quadratic Zk,2 parts of Z^ satisfy the recursion systems 

Zk+i,i = dZk,i, Z^,^ = U, (4.28) 

fc-i 

Zk+1,2 = dZk,2 + ^ Zm,lZk-m,l, ^3,2 = t^2- (4.29) 

m=l 

Mathematical induction gives (4.27). 
Conclusion 4.4. There holds for k > 2 

Gk = U2k-2 + { ('^^^ + (-1)-} U,U2k-.-p + ('^^) U\ (4.30) 
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where the points indicate terms of degrees higher than two and less than k. 

Proof. The linear i and the highest degree Gk,k parts of Gi~ follow immediately from 

(4.11). Therefrom we have also for the quadratic part Gk,2 

fc-i 

Gk,2 = 2 ^ Z2m-l,lGk-m,l + -^2fc-l,2- (4-31) 
m=l 

Explicit calculation gives the result. 

Note that in the scalar case Gk,2 had been found by means of a completely different 
method by Gilkey^^. 

Let us now derive another recursion system for the differential polynomials G^ ■ Let in 
the following the square and the curly brackets denote the commutator and anticommutator 
of matrices respectively. 

Theorem 4.5. There exist unique differential polynomials 

Wk=:d-^[Gk,U], (4.32) 

to the sequence Gk{k = 1,2, . . .) of differential polynomials in U defined by (4.11), (4.8), 
(4.9), so that 

[Gk,U] = dWk, (4.33) 
and the G^ satisfy the differential-recursion system 

dGk+i = d^Gk + {Ud + dU,Gk} + [d-^[Gk,U],U]. (4.34) 

Proof. Let us consider once more the equation for the resolvent kernel 

d^G{x, x'; A) + {U{x) - XI)G{x, x'; A) = -6{x - x')I. (4.35) 

As assumed from the beginning, the matrix-valued potential is Hermitian, = t/, so that 
the Schrodinger operator is self-adjoint, — L, and the resolvent satisfies the symmetry 
condition 

g'^{x,x';X) = G{x',x;X). (4.36) 

Using this property and interchanging the variables x and x' we get another differential 
equation for G: 

d'^G{x, x'; A) + G{x, x'; X){U{x) - XI) = -d{x - x')I. (4.37) 

where d' = d/dx' . 

The difference of (4.35) and (4.37) reads 

{d + d'){d - d')Gix, x'; A) = Gix, x'; X)Uix') - U{x)G{x, x'; A). (4.38) 

For every smooth two-point function H{x, x') there holds 

[{d + d')H] = d[H]. (4.39) 
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Using this rule and restricting (4.38) to the diagonal x = x', we get 

[[G],U]^d[{d-d')G], (4.40) 

or put in another way 

d-^[[G],U] = [{d-d')G]. (4.41) 

Both [G'(A)] and W{\) := [{d — d')G{X)] admit asymptotic expansions in powers of A"-*^/^, 
where the coefRcients Gk and Wt respectively are differential polynomials in U. Hence 
there are unique differential polynomials Wk such that 

dWk = [Gk,U]. (4.42) 

Now, let us apply the operators {d + 3d') and {d' + 3d) to the equations (4.35) and 
(4.37) respectively and then take the sum. Restricting to the diagonal the result 

{d + d')^G{x, x'; A) - 4X{d + d')G{x, x'; A) 

+ {dU{x))G{x, x'; A) + G{x, x'; \)d'U{x') 

+ 2U{x) {d + d')G{x, x; A) + 2{{d + d')G{x, x'; X))U{x') 

- U{x){d - d')G{x,x'; A) + {{d - d')G{x,x'; X))U{x') = 0, (4.43) 

and using (4.39), we obtain 

d^[G] - 4Xd[G] + 2{U, d[G]} + {dU, [G]} + [d-^[[G],U],U] = (4.44) 

Substitution of the asymptotic expansion of the resolvent kernel on the diagonal (4.25) 
proves that the quantities Gk satisfy the recursion system (4.34). 

The formula (4.34) has been found first by Olmedilla et al.^°, (see also the recent 
paper of Bilal^-*^). In the scalar case, when all commutators vanish, it is known under the 
name 'Lenard recursion'. As a matter of fact, Burchnall and Chaundy^'^ discovered the 
sequence {Gk) and the scalar 'Lenard recursion' already in the twenties. 

Let us write (4.34) in the form 

dGk+i = BGk (4.45) 
where the linear operator B is defined by 

EH := d^H + 2{U, dH] + {dU, H} + [d'^ [H, U] , U] (4.46) 

Since the operator d is formally anti-selfadjoint, d^ — —d, so is B: 

st = -B. (4.47) 
In the scalar case all the commutators vanish and B becomes a differential operator 

B = d^ + 2Ud + 2dU, (4.48) 
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where the factor U means multiphcation by U. The recursion system (4.45) can be formally 
solved: 

Gk = l{d-'B)'^-I. (4.49) 

Thereby, the action of the positive powers of the operator B on the identity matrix is 
well-defined and produces the differential polynomials Gk- 

A special motivation for the study of the one-dimensional Schrodinger equation is 
its relation to the Korteweg-de Vries (KdV) hierarchy. If one assumes a one-parameter 
dependence of the potential U = U{x,t), then the differential polynomials — Gk[U]{k = 
1,2,...) constructed above constitute the KdV hierarchy in the scalar case: 

^ ^^G,[U]. (4.50) 

For k = 1 a, linear differential equation emerges, for k = 2 we have the KdV equation itself, 
for k > 2 the equation (4.50) is called a higher KdV equation. Moreover, it is known that 
(in the scalar case again) the quantities Gk as well as the are conserved densities to 
the KdV equation. Several papers discuss a matrix KdV equation^^'^^'^^ 

—U = —G2 = {d^fU + Z{d,U)U + W{d,U). (4.51) 

OT OX 

In particular, Olmedilla^^ proved that (4.51) is solitonic in the sense that the inverse 
scattering method works for (4.51). One might guess a solitonic character for the whole 
matrix KdV hierarchy (4.50) (/c = 2, 3, . . .). 
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